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Abstract. The three— dimensional incompressible viscous Boussinesq equations, under the assumption of hydro- 
static balance, govern the large scale dynamics of atmospheric and oceanic motion, and are commonly called 
the primitive equations. To overcome the turbulence mixing a partial vertical diffusion is usually added to the 
temperature advection (or density stratification) equation. In this paper we prove the global regularity of strong 
solutions to this model in a three-dimensional infinite horizontal channel, subject to periodic boundary conditions 
in the horizontal directions, and with no-penetration and stress-free boundary conditions on the solid, top and 
bottom, boundaries. Specifically, we show that short time strong solutions to the above problem exist globally in 
time, and that they depend continuously on the initial data. 
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1. Introduction 

The partial differential equation model that describes convective flow in ocean dynamics is known to be the 
Boussinesq equations, which are the Navier-Stokes equations (NSE) of incompressible flows with rotation coupled 
to the heat (or density stratification) and salinity transport equations. The questions of the global well-posedness 
of the 31? Navier-Stokes equations are considered to be among the most challenging mathematical problems. In 
the context of the atmosphere and the ocean circulation dynamics geophysicists take advantage of the shallowness 
of the oceans and the atmosphere to simplify the Boussinesq equations by modeling the vertical motion with the 
hydrostatic balance. This leads to the well-known primitive equations for ocean and atmosphere dynamics (see, 
e.g., [24], [25], [29], [31], [32], [37], [38] and references therein). A vertical heat diffusivity is usually added as a 
leading order approximation to the effect of micro-scale turbulence mixing (cf., e.g., [16], [17], [24], [31]). As a 
result one arrives to the following dimensionless 31? variant of the primitive equations (Boussinesq equations): 

dv Qv -» 

— + (v-V H )v + w— + f o kxv + V H p + Liv = (1) 
at oz 

d zP + T = 0, (2) 
V H ■ v + d z w = 0, (3) 
dT dT 

— +vV H T + w—+L 2 T = Q, (4) 
at oz 

where the horizontal velocity vector field v = (vi, v 2 ), the velocity vector field (v%, v 2 ,w), the temperature T and 
the pressure p are the unknowns. f is the Coriolis parameter, Q is a given heat source. For simplicity, we drop 
the coupling with the salinity equation, which is an advection diffusion equation, but the results reported here 
will be equally valid with the addition of the coupling with the salinity. Moreover, we also assume for simplicity 
that Q is time independent. The viscosity and the heat vertical diffusion operators L\ and L 2 , respectively, are 
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given by 

1 a 1 ° 2 

L ' = -R^-R^ (5) 
1 d 2 

L2 = 'If, W (6) 
where J?i,J?2 are positive constants representing the horizontal and vertical dimensionless Reynolds numbers, 
respectively, and R3 is positive constant which stands for the vertical dimensionless eddy heat diffusivity turbu- 
lence mixing coefficient (cf., e.g., [16], [17]). We set Vjj = (d x ,d y ) to be the horizontal gradient operator and 
Ah = d 2 + d 2 to be the horizontal Laplacian. We denote by 

T u = {(x,y,0) el 3 }, (7) 

T b = {(x,y,-h)eR 3 }, (8) 

the upper and lower solid boundaries, respectively. We equip system (1)— (4), on the physical top and bottom 
boundaries, with the following no-normal flow and stress free boundary conditions for the flow velocity vector 
field (v,w), namely, 

dv 

on Y u : — = 0, w = 0, (9) 
dv 

on T fc : — = 0, w = 0, (10) 

and for simplicity, we set the Dirichlet boundary condition for T: 

71 Z=0 = 0, T\ z= _ h = l. (11) 
Horizontally, we set (v, w) and T to satisfy periodic boundary conditions: 

v(x + 1, y, z) = v(x, y + 1, z) = v(x, y, z)\ (12) 
w(x + 1, y, z) = w(x, y + l,z) = w(x, y, z); (13) 
T(x + 1, y, z) = T(x, y + l,z)= T(x, y, z). (14) 

We will denote by 

M = (0, l) 2 and O = Mx(-/i,0). 
In addition, we supply the system with the initial condition: 

v(x,y,z,0) = v {x,y,z), (15) 

T{x,y,z,0) = T {x,y,z). (16) 

System (1)— (16) is a modified form of the rotational Rayleigh-Benard convection problem taking into consideration 
the geophysical situation of the shallowness of oceans and atmosphere. The original three-dimensional Rayleigh- 
Benard convection model (which is identical, in the absence of heat diffusion, to the Boussinesq model of stratified 
fluid) has been a subject to study for many years, numerically, experimentally and analytically (see, e.g., [3], [8], 
[12], [18], [28], [29], [30], [37], and references therein). However, the question of global regularity is still open and is 
as challenging as the 31? NSE. Recently, the authors of [9] and [20] have shown the global well-posedness to the 2D 
Boussinesq equations without diffusivity in the heat transport equation (see also, [11], for recent improvement). 
In [7] it is observed that thanks to hydrostatic balance (2) the unknown pressure is essentially a function of the 
two-dimensional horizontal variables. We take advantage of this observation in [7] to establish the L 6 estimates 
for the velocity vector field which allows us to prove the global well-posedness of the 3D primitive equations under 
the geophysical boundary conditions. In [22] the authors take advantage of this observation as well, and proved 
the global well-posedness of the 3D primitive equations with the Dirichlet boundary conditions by dealing directly 
with the "pressure" which is a function of two variables. In this paper we study system (1)— (16), exploring again 
the hydrostatic balance which leads to an unknown "pressure" that is a function of only two variables, and use 
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the techniques and ideas introduced in [7], [9] and [20], to show in section 3 that strong solutions exist globally 
in time provided they exist for a short interval of time. Furthermore, we show in section 4 the uniqueness and 
continuous dependence on initial data of these strong solutions. The short time existence of strong solutions to 
this model will be reported in a forthcoming paper. 

This paper is organized as follows. In section 2, we reformulate system (1)— (16) and introduce our notations 
and recall some well-known inequalities. Section 3 is the main section in which we establish the required estimates 
for proving the global existence in time for any initial data. In section 4 we prove the uniqueness of the solutions 
and their continuous dependence on initial data. 

2. Functional setting and Formulation 
2.1. Equivalent Formulation. We denote by 

i r° 

<t>{x, y) = - J cj>{x, y, z)dz, V (x, y) G M; (17) 

and denote the fluctuation by 

= 0-£ (18) 

Notice that 

4> = Q. (19) 
Similar to [7], by integrating (2) and (3) vertically, we get 



w 

and 



(x,y,z,t) = - [ V H -v(x,y,Z,t)dt, (20) 
J-h 



p(x,y,z,t) = - [ T(x,y,£,t)d4 +p a (x,y,t), 

J-h 



(21) 



where p s is the pressure on the bottom z = —h. Essentially, p s (x,y,t) is the unknown pressure, and we observe, 
as before, that it is a function of two spatial variables (x, y). As we mentioned in the introduction we explore this 
property as in [7] (see also [22]) to prove our global regularity result. 

Replacing T by T + jfc, we have the following equivalent formulation for system (1)— (16): 

— +L 1 v + (v V H )v - ( y V ff ■ v(x, y, £, i)d£ J — 



h 



+f kxv + V H Ps(x,y,t)-V H I T(x,y,€,t)d£ = 0, (22) 

J-h 

V H -v = 0, (23) 

^+L 2 T + v V H T - ^J" \7 H ■ v(x, y, C, t)<%) + ^) = <3. (24) 



dv 
~d~z 



dv 



z=0 9Z 



= 0, v(x + l,y,z) = v(x,y + l,z)=v(x,y,z), (25) 

Z — — h 

T\ z=0 =T\ z= _ h = 0, T(x+l,y,z) = T(x,y + l,z)=T(x,y,z), (26) 

v(x,y,z,0) = v (x,y,z), (27) 

T(x,y,z,0) = T (x,y,z)-^. (28) 
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In addition, v and v satisfy the following coupled system of equations: 
dv 1 



— - -j^^HV + (v ■ V H )v + [{v ■ V H )v + (Vjj • v) v\ 



+f k xv + \7 H 



p s (x,y,t)-- / T(x,y,£,t) d£ dz 
n J~hJ-h 



(29) 



-hJ-h 

V H -v = 0, (30) 
-^+L 1 v + {v-V H )v- yj V H -v(x,y,£,t)d£J + (v • V H )v + (v ■ V H )v 

-[(v-V h )v + (Vh-v)^ + f k x v - V H [J ^ T(ar, y, £, t)d£ ~\ J J T ^ V> & = °- ( 31 ) 

2.2. Functional spaces and inequalities. Let us denote by L"(fl),L q (M) and W m ' q (fl), W m ' q (M), and 
H m (Q) =: VF m ' 2 (fi),iJ m (M) =: VK m < 2 (M), the usual L«-Lebesgue and Sobolev spaces, respectively ([1]). We 
denote by 

||0|| =/ (JpJ^^)! 9 dxdydz)~ q , for every e L 9 (ft) ^ 

' \ (/ M |^(a:,i/)|«da:dy)«, for every e L«(M). 

For convenience, we recall the following Sobolev and Ladyzhenskaya type inequalities in (IR/Z) 2 and in f2 (see, 
e.g., [1], [10], [15], [23]) 

U\\l* { m) <C \\<l>\)$\\<l>t£ m , tyeH^M), (33) 

II0I|ls ( m) <C ||<A||^ M) ||<A||^ 4 (M) , Vtfei^M), (34) 

||V ff 0|| L 4 (M) <Co||0||L /2 ||^ll^ 2 (M)' V^ei/ 2 (M), (35) 

||V ff 0|| L 4 (M) < C o ||0||^ 2 (M) ||V ff 0||L /2 + IHI^m), V <$> such that V ff e L°°(M), (36) 



and 



HV-ll^nj < CollVll^nj HV-ll^n), (37) 



U«(n) <Co||^|| H i(o), (38) 

for every V € i? 1 (fi). Here Co is a positive scale invariant constant. Also, we recall the following version of 
Helmholtz-Weyl decomposition Theorem (cf. for example, [26], [15], [39]) 

||V^|k-,, (M) < C (||V ff • <j>\\ W m, q{M) + ||V£ ■ cj>\\ W m, HM) ) , (39) 



2 - 

for every <j> G (W m ' q (Mj) . Moreover, we recall the following Brezis-Gallouet or, Brezis-Wainger inequality (see, 
e.g., [26], [4], [5], [13]) ' 

||<AIU~(M) < CU\\ HHM) (l + log+||0|| ff 2 (M) ) 1/2 , (40) 

for every <j> G H 2 (M), where log + r = logr, when r > 1, and log + r = 0, when r < 1. Also, we recall the following 
inequality (see, e.g., [2] and [21]) 

||V H 0|| L oc (M) < C(||V H • <^|| L oc (M) + || V ff x <p\\ L oo {M) ) (1 + log+ || Vh4>\\hhm)) , (41) 

for every V H 4> G H 2 (M). Moreover, by (33) we get 

\\K\m) = II W \\i H M) < C\\ \4>\ q \\h m \\ W \\ 2 m( M) 

< CMlll (/ m l^l 29 " 2 |V«0| 2 dxdy^j + U\\tl (42) 
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for every (f> satisfying J M \<f>\ 2q ~ 2 |V#0| 2 dxdy < oo and q > 1. Also, we recall the integral version of Minkowsky 
inequality for the L p spaces, p > 1. Let Sli C W ni and f2 2 C l™ 2 be two measurable sets, where mi and m 2 arc 
two positive integers. Suppose that f(£,r]) is a measurable function over f!i x O2. Then, 



< 



1/p 



(43) 



Finally, we recall the following inequality from Proposition 2.2 in [6] 



/(/" 

JM \J-h 



ip!(x,y,z) dz 



ip 2 (x,y,z)ip 3 (x,y,z) dz) dxdy 



< C||Vi||^ /2 ||VhVi||2 /2 IIV' 2 ||2 /2 I|VhV2||2 /2 ||V'3||2 + IIV1II2IIV2II2IIV3II2, 

for every V^V^ G i? 1 ^) and ^3 G 7 2 (^), and 

J i^J ipi(x,y,z) dz^j (^j \\7 H ip2(x,y,z)\ip 3 (x,y,z) dz^j dxdy 

< C||^i||2 /2 ||V^i|| 2 /2 ||^||^ 2 ||V ff V^ 2 ||^ /2 ||^3||2 + HV1II2IIV2II2IIV3II2, 
for every Vi G tf 1 ^), V H ^ 2 € ff 1 ^) and ^3 G 7 2 (^)- 

3. Global existence of strong solutions 



(44) 



(45) 



In the previous section we have reformulated system (1)— (16) to be equivalent to system (22)-(28). In this 
section we will show that strong solutions to system (22)-(28) exist globally in time provided they exist in short 
time intervals. 

Theorem 1. Let Q G H 2 (n),vo G H 4 (n),To G H 2 {Vl) and T > 0. Suppose that there exists a strong solution 
(v(t),T(t)) of system (22)-(28) on [0,7] corresponding to the initial data (vo,To) such that 

A H v z , V H TeC([0, T},H\n)), 

v zz ,A H V H v z , V H T Z G L 2 {[0,T\,H l {Q)). 

Then this strong solution (v(t),T(t)) exists globally in time. 

Remark 1. Notice that one can recover the pressure p s from system (29)-(30) in the same way as in 2D NSE 
(see, e.g., [10], [34], [35]). 

Proof. Let [0,%) be the maximal interval of existence of a strong solution (v(t),T(t)). In order to establish the 
global existence, we need to show That % = 00. If % < 00 we will show ||Ajji> z (i)||jyi(£2), ||VfjT(£)||#i(o), 
ds,J>\\V H T z (s)\\ 

ff!(fi) an d Jo H t ' Z2 ( s )llff 1 (n) are all bounded uniformly in time, for 
i G [0,T*). As a result the interval [0,%) can not be a maximal interval of existence, and consequently the strong 
solution (v(t),T(t)) exists globally in time. 

Therefore, we focus our discussion below on the interval [0,71). 

3.1. Hull! + \\TW2 estimates. By taking the inner product of equation (24) with T, in L 2 (Q), we get 

ldjjTJH 1 2 
2 dt + D Uzh 



R 3 



= QT dxdydz - j v • V ff T - V ff • v(x, y, £, t)d^j + \ 



T dxdydz. 



Integrating by parts and using the boundary condition (26), we get 



~ L { v ' VhT ~ (/ h Vh ' y ' ^ * )d © T dxdydz = 



(46) 



(47) 
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As a result of the above we conclude that 
1 



2 dt 



R 3 



mi 



T dxdydz < \\Q\\ 2 ||T|| 2 + || V H «|| 2 ||T|| 2 . 



(48) 



Moreover, by taking the inner product of equation (22) with v, in L 2 (il), we reach 

MMIi , i ..- - i 



2 



|Vffv||? + ^-K||i 



= - J (v- V H )v - (^j V H ■ v{x, y, £, t)d£j |^ 



v dxdydz 



In ( 



w dxdydz 



First, we notice that 



(49) 



(50) 



Next, by integration by parts and using the boundary conditions (25), in particular, the horizontal periodic 
boundary conditions, we get 

dv 



dz 



v dxdydz = 0. 



(51) 



/ (v-Vh)v-([ V H -v(x,y,Z,t)dt 
Jn L \J-h 

Thanks to (30) and, again, the horizontal periodic boundary conditions, we also have 

/ ^HPs{x,y,t) ■ v(x,y,z,t) dxdydz — h / V hPs -v dxdy = -h I p s (V H -v) dxdy = 0. (52) 
Jn Jm Jn 



in Jm Jn 

By integration by parts, the periodic boundary conditions (25), and Cauchy-Schwarz inequality, we obtain 



' H 



T(x, y, £, t)d£ I • v dxdydz 



<h\\T\\ 2 \\V H v\\ 2 . 



(53) 



Thus, by (48)-(53) we have 

MINi + imil) . 1 

2 



dt 



+ ^Hv\\l + ^\v z \\i + ^-\\T z \\l 



< ||Q|| 2 ||T|| 2 + (l + /i)||T|| 2 ||V ff w|| 2 . 
By Cauchy-Schwarz inequality, we obtain 

d(\\v\\ 2 2 + \\T\\ 2 2 ) 1 2 1 2 1 2 

dt " 2 " 2 S^" " 2 

<||Q||i + (l + i2i)(l + /i) 2 ||r||l. 
Thanks to Gronwall's inequality we get 

\Ht)\\l + \\T(t)\\l < C {\\voWl + \\T \\ 2 2 ) ed+^H^) 2 * + C\\Q\\ 2 2 ; 



and 



/' 

Jo 



Ri 



\Vhv{s)\\1 + —\\v x ( 



+ 



\Tz{a) 



ds<C 



(lMl + ||T ||l) 



WQWh 



(54) 

(55) 
(56) 

(57) 
(58) 
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Therefore, for every t £ [0,%), we have 

IK*)II! + + f \\\Vhv{s)\\1 + \\v z (s)\\l + \\T z (s)\\l] ds < K U (59) 

Jo 



where 



K X = C 



IkolU + llTolU) e (l+^)d+^) 2 * + HQIII J . (60) 



3.2. 1 1 r || oo estimates. We follow here the idea of Stampaccia for proving the Maximum Principle. The proof we 
present here is also similar to the one in [14] (see also [36]). Denote by r(t) = T{t) — (1 + ||Xo||co + ||Q||oo t). It is 
clear that t satisfies: 

-^+v-\7 h t + w— + L 2 t = Q-\\Q\\ co . (61) 

Let t + = max{0, r} which belongs to and satisfies 

t + (z = 0) =t + {z= -h) =0. (62) 

Taking the inner product of the equation (61) with r + in L 2 (f2) and applying the boundary conditions (62), we 
get 

^<3l3l + 1_ \\d z T+\\l = [ {Q - ||Q||oo)t + dxdydz - [ [v ■ V h t + w 8 z t] t+ dxdydz. (63) 
By integration by parts and using the boundary conditions (26) and (62), we get 

/ [v ■ V h t + w 8 z t] t+ dxdydz = 0. (64) 

JQ. 



Thus, 



ld\\r+\\l , 1 



+ ^\\d z T+\\l dxdydz = [ (Q- ||Q||oc)t+ dxdydz < 0. (65) 

^3 JO 



2 dt R 3 
Therefore, we obtain 

\\r + (t)\\t<\\r + (t = 0)||1 = 0. (66) 

Thus, r + (t) = 0. As a result, we have 

T(t)<l+||2b||oo + ||Q||oot- (67) 
By applying similar arguments, we also have 

r(t) > -(1 + ||3b||oo + IIQIloo *)■ (68) 
Therefore, T satisfies the following estimate: 

ITOIIoc < K 2 = 1 + IITolU + IIQHoo t. (69) 
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3.3. 1 1 v 1 1 6 estimates. Taking the inner product of the equation (31) with \v\ 4 v in L 2 (tt) and using the boundary 
conditions (25), we get 

IdNIi , 1 



6 dt 



\v z \ 2 \v\ 4 + d z \v\ 2 r\v\ 2 ) dxdydz 



+ ~kJn ( |V ^ 21 ^ 4 + I Vff I"' 2 ! 2 l^l 2 ) dxd V dz + J (l 



+/ fcx«-V H ( / T(x,y,£,t)d£-- 



hJ-h 



T(x,y,£,t)d£dz I > • \v\ 4 v dxdydz. 



Observe, again, that 

(/ofc x v) • |w| 4 w = 0. 
Moreover, by integration by parts and the boundary conditions (25), we also get 



/ (v-Vh)v-I V H -v(x,y,£, 
Jn l \J-h 



dv 



dz 



\v\ 4 v dxdydz = 0. 



(70) 



(71) 



Furthermore, by virtue of (30) and by the boundary conditions (25), in particular the horizontal periodic boundary 
conditions, we have 



/ (v(x,y,t) ■ \7 H )v(x,y,z,t) ■ \v(x, y, z, t)\ 4 v(x, y, z, t) dxdydz = 0. 
Jn 



(72) 



Thus, (70)-(72) imply 
6 dt Ri 



l n (\V H v\ 2 \v\ A +\V H \v\ 2 \ 2 m 2 ) dxdydz +±- 



v z \ 2 \v\ 4 + 



\d z \v\ 2 \ 2 \v\ 2 ^j dxdydz 



= - J | {v ■ V H )v - (v ■ Vh)v + (Vh • v) v 



-Vh ( / T(x,y,£,t)dS-- 



hJ-h 



T(x,y,^,t)d^dz > I • \v\ 4 v dxdydz. 



(73) 



Notice that by integration by parts and using the boundary conditions (25), in particular the horizontal periodic 
boundary conditions, we have 



- / (v ■ V H )v - [(v ■ Vh)v + (Vh • v) v\ 
Jn 1 



(74) 



r z 



-V fl ( / T{x,y,Z,t)d£-- 



T(x,y,£,t)d£dz 



hJ-h 



\v\ 4 v dxdydz 



2 3 



(V H • v) v ■ \v\ 4 v + (v ■ V H )(\v\ 4 v) ■ v - J2J2v k v j d Xk {\v\ 4 v j ) 

fc=i j=i 

r z 



T(x,y,S,t)dt-l f f T(x,y,^t)d^dz)v H -(\v\ 4 v 
h n J-hJ-h / 



T) 



dxdydz. 



(75) 
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As a result, we obtain 

IIS 



ld\\v\\ e 



6 dt Ri 

< C 



kL ( |V ^ |2| ^ |4+ l V ^| 2 | 2 |^ 2 ) dxdydz+-^J^(\v z \ 2 \v\ 4 +\d z \v\ 2 \ 2 \v\ 2 ) dxdydz 



M 







\\7 H v\ \v\ 5 dz 



dxdy + cj ^ \v\ 2 dz)(^J \V H v\\v 



~i i~|4 



dz 



dxdy 



+C ( \\T\ f \V H v\ \v\ a dz 

JM L J-h 



dxdy. 



Therefore, by the Cauchy-Schwarz inequality and Holder inequality we reach 



ld\\v\\\ 



6 dt 

< 



-^-J (\V H v\ 2 \v\ 4 + \S7h\v\ 2 \ 2 \v\ 2 ) dxdydz +^-J^{\v z \ 2 \v\ 4 +\d z \v\ 2 \ 2 \v\ 2 ) dxdydz 



C f \v\( f \Vhv\ 2 |«| 4 dz) ' ( f \v\ 6 dz 

JM [ \J-h / \J-h 

+C J \v\ 2 dz) (^J° \\7 H v\ 2 \v\ 4 dz) ' [J \vt dz 



dxdy 



1/2 



dxdy 



+C 



A I 



\T\ ( / \V H v\ 2 \v\ 4 dz 

h 



1/2 , -0 



u| 4 dz 



1/2 



dxdy 



Moreover, 



6 ' i?i 

< C\\v\\ L i {M) 



kL ( |V ^ |2|?|4+ l Vff l^ 2 | 2 ^ 2 ) dxdydz + J n (l^| 2 |^| 4 + \ d M 2 \ 2 \v\ 2 ) dxdydz 



\V H v\ 2 \v\ 4 dxdydz) ' (j (^J \v\ G dz) dxdy) 



+C 



M 



1/4 



v\ dz] dxdy 



a 



|V ff u| 2 |u| 4 dxdydz 



|T| || L 4 (M) (jf |V H «| 2 |«| 4 dxdyd^ 7 ^ ^ |«| 4 dz) dxdy 
Using the Minkowsky inequality (43) we get 

( I ( I Fl 6 dz ) dxd v) < C j ( [ \v\ 12 dxdy 

\JM \J-h ) j J-h \JM 



1/2 



M 





-h 
1/4 



1/4 



v\ dz] dxdy 



1/2 



dz. 



Thanks to (42), 



M 



\v\ 12 dxdy < Co ( f \v\ 6 dxdy) ( f \v\ 4 \V H v\ 2 dxdy) + ( f \v\ 6 dxdy 

\Jm J \Jm J \Jm 



Thus, by the Cauchy-Schwarz inequality we obtain 

2 \ 1/2 



^ \vf dz) dxdyj < C\\v\\ 3 L e {n) \v\ 4 \V H v\ 2 dxdydt 



1/2 



+ v 



L6(f2)- 
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Similarly, by (43) and (34), we also obtain 

2 \ 1/2 



and 



(/ (/ !^ 4dz ) dxdy j - C J (/ fi\ 8 dxd v) dz 

< C [ \\v\\h m (\\V H v\\ LHM) + \\v\\ LHM) ) dz < C\Hl (\\V H v\\ 2 + \\v\\ 2 ) . 

J —h 

1(1 \v\ 2 dz] dxdy] < C [ ( f |w| 8 dxdy] ' dz 

JM \J-h J J J-h \JM J 

f M%\m) (||VHI^(m) + ^LHM)) dz < C\\vf 6 /2 (\\V H v\\ 2 + Wvh) 1 

J — h 



(82) 



< C 



Therefore, using (81)-(83) and (33), we reach to 



(83) 



1« , 

6 dt Ri 

< C\\v\\\' 2 \\V H v 



^ J^(\V h v\ 2 \v\ 4 + \Vh\v\ 2 \ 2 \v\ 2 ) dxdydz + ^-J^{\v z \ 2 \v\ 4 +\d z \v\ 2 \ 2 \v\ 2 ) dxdydz 



II 1/2 1. — ..3/2 



\V H v\ 2 \v\ 4 dxdydz 



3/4 



C\\v\\l{\\V H v\\ 



+C\\v\\\ /2 llVHI^II^IIe + CIITl^ \\v\\l /2 {\\V H v\\ 



\ ,2 +\\v 



w /2 ) 



IN.)(/ n 

/ |V#w| 2 \v\ 4 dxdydz 
Jn 



\S7 H v\ 2 \v\ 4 dxdydz 



1/2 



1/2 



By Young's inequality and Cauchy-Schwarz inequality we have 



^ ' I (\^HV\ 2 \V\ 4 



dt 



Ri 



V H \v\ 2 \ 2 \v\ 2 ) dxdydz + —j (\v z \ 2 \v\ 4 +\d z \v\ 2 \ 2 \v\ 2 ^j dxdydz 



< C\\v\\l \\V H v\\l\\v\\ 6 e + C\\v\\l\\V H v\\l + CWT^ + C\\v\\ 2 \\v\\l 
Thanks to (59), (69) and Gronwall inequality, we get 

\\v(t)f 6 + f f-l- / \V H v\ 2 \v\ 4 dxdydz + 1- f \v z \ 2 \v z \ 4 dxdydz) < K 3 , 
Jo \ H i Jn K 2 Jn J 



where 



K 3 = e 



Kit 



lkl| 6 ffl( n)+^2 4 i 



(84) 



(85) 



3.4. || Vjru||2 estimates. By taking the inner product of equation (29) with -A^jj in L 2 (M), and applying (30), 
and using the boundary conditions (25), we reach 

ld||V ff iJ||| , 1 



dt 



Ri 



\\^HV\\ 2 2 



(v ■ Vh)v + [(v ■ V h)v + (Vh ■ v)v] + f k x v j • A H v dxdy. 



Following the situation for the 2D Navier-Stokes equations (cf. e.g., [10], [32]) we have 



/( 

Jm 



(v ■ Vij)» • Ajjv dxdy 



< C||iJ||^ /2 ||V ff U|| 2 \\A H v 



1 1 3/2 



(86) 



(87) 
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By the Cauchy-Schwarz and the Holder inequalities, we have 



/ (v ■ Vjj)d + (Vij • v) v ■ Ahv dxdy < C / \v\ \Vhv\ dz \ Ahv\ dxdy 
Jm Jm J-h 

cj (^j ^\v\ 2 \S7 H v\dz^j (^j \V H v\dz^j \A H v 



< 



< C 



M 



< C||V^||^ /2 



\v\ \V hv\ dz I dxdy 



\v\ 4 \V H v\ 2 dxdydz 





1/4 






1/4 




dxdy 




Im ( 


/ |Vjyt)| dz I dxdy 




/ \A H v\ 2 dxdy 












Jm 



1/2 



1/4 



\A„v 



HV\\2, 



and 



/ foxv- A H v dxdy 
Jm 



<C\\v\\ 2 \\A H v\\ 2 . 



Thus, by Young's inequality and the Cauchy-Schwarz inequality, we have 

^JfM + ±-\\A H V\\l < C\\v\\l\\V H m + C\\V H v\\ 2 2 +C [ \v\ 4 \V H v\ 2 dxdydz + C\\v 
at Hi j n 

By (59), (84) and thanks to Gronwall inequality we obtain 

\\Vhv\\ 2 2 + 4- f \A H v\l ds < K A , 



-l|2 
2- 



where 



K 4 = e 



Kit 



\vo\\m ( n)+K2 + K 3 



(88) 
(89) 
(90) 

(91) 
(92) 



3.5. \\v z \\e estimates. The a priori estimates (59)-(91) are essentially similar to those obtained in [7]. From now 
on, we will get new a priori estimates of various norms. 
Denote by u = v z . It is clear that u satisfies 



du 

— + Liu + (v ■ V H )u 



-(u ■ V H )v - (V ff • v)u + f k x u - V H T = 0. 



= -h 



0. 



(93) 
(94) 



Taking the inner product of the equation (93) with m|u| 4 in L 2 , we get 



6 dt Ri 



R 2 ' 



Notice again that 



= I i} V ' ^ H ^ U ~ (/ Vh ' V<yX ' Vl ^ Ijz) ' "' U ' 4 dxdydz 
- J ({u ■ Vh)v - (Vh • v)u + fok x u- VrT^ ■ u\u\ 4 dxdydz. 

fok x u ■ u|w| 4 = 0. 



(95) 
(96) 
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Integrating by parts and using the boundary conditions, in particular (94), give 

du 



(v ■ V H )u 



V h ■ v(x,y,£,t)d£ 



dz 



u\u\ dxdydz = 0. 



(97) 



Thus, by (96), (97) and Holder inequality, we have 



idH| 

6 dt 



■jrW h 2 \v H u\ h» + A|| H 2 ldzU \ ni 

til X12 



< 



— / ((u ■ V h)v — (Vh • u)it — V ijT) • u|u| 4 dxdydz 
Jn 

C [ \v\ \u\ 5 \V H u\ dxdydz + C [ \T\ |u| 4 |V#it| ctedycte 
in in 



<C(|M| 6 || U 3 || 3 ||V ffU 2 || 2 + !|T|| 6 \\u%\\V H u%) 
< C (\\v\\ 6 \\u\\l /2 \\V H u 3 f 2 /2 + \\T\\ 6 \\u\\l\\V H u%) . 
Thanks to the Cauchy-Schwarz inequality, we have 



(98) 



d\\u\\ 
dt 



s-W H 2 l v ^l II2 + -5-II M 2 \d*u 



Ri 



R2 



<C(l + \\v\\t) ||«||g + ||T||{ 
<C(l + \\\7 H vf 2 + \\v\ 



\\< + \\nt 



Using (59), (84), (91), and Gronwall inequality, we get 

rt r 



Il«ll2 + 



/ 

Jo 



^11 H 2 Whu\ + M 2 \d z u\ \\l 
tx\ til 



ds < K 5 , 



where 



\\dM\ & H i m +K & 2 t 



(99) 
(100) 

(101) 
(102) 



3.6. ||w zz ||2 estimates. Taking the inner product of the equation (93) with — u zz in L 2 {Q) and recalling that 
u = v z , which satisfies the boundary condition (94), we get 



/ II 2 



dt 



■-E-\\Vhu x \\1 + -±-\\u xx \\1 



du\ 



V ff • v(x, y, £, i)d£ ) —j ■ u zz dxdydz 



= 1 ((vvh)u- ( r 

Jn \ \J-h 
+ J {{u ■ V//)u - (Vh • v)u + fok x u — VrT^J ■ u zz dxdydz 

= - [(u- Vh)u + (u z ■ Vh)v + (u ■ Vh)u - (Vh • u)u - (Vh • v)u z ] ■ u z dxdydz - j 
Jn Jn 

< CH| 6 ||Vhu|| 2 |K||3 + C||w||6||VffU z || 2 ||«z||3 + ||T,||2||Vir« z ||2 

' +||T z || 2 ||V ff «,||2. 



T z (Vh ■ u z ) dxdydz 



,1/2 



<C[||u|| 6 ||Vffu||2 + ||«|| 6 ||Vff« 2 || 2 ]K|| 2 
By the Cauchy-Schwarz and Young's inequalities, we have 

id|M| 1 2 1 



V H u z \\\ /2 + \\u zz \\\ /2 



< C (\\v\\i + \\u\\i) \\u z g + C\\V H u\\ 2 + C\\T z \\l 
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Applying (59), (84), (101), and Gronwall inequality yield 



f 

Jo 



Ri 



R.2 



ds < Kq, 



where 



K 6 = Ce^3 /3 +^ 2/3 )* 



\vo\\h(n)+ K i 



(103) 



(104) 



3.7. ||V# x v z \\% + ||Vif • v z + RiT\\l estimates. Let /? be the solution of the following two-dimensional elliptic 
problem with periodic boundary conditions: 

A H /3 = V H T, [ I3dxdy = 0, (105) 

where z is considered as a parameter. Roughly speaking, /3 is like the potential vorticity. Notice that 

V ff • p = T, V ff x/J = 0. (106) 

Recall that u — v z . We denote by 

( = u + R 1 f3, (107) 

V= (Vir -C) = -u = d x u 2 -d y u u (108) 
= (V H ■ C) = V ff • u + RiT = d xUl + d y u 2 + RiT. (109) 
Applying (39) for ( and (i with m > 0, 1 < q < oo, using (106)-(109), we obtain 

||Vjju||^m,«(M) < C (\\V H (\\w™>«(M) + Ri\\VhP\\w™-«(m)) 

< C (\\v\\w^HM) + PWw^HM) + \\T\\ wm . HM) ) , (110) 
where, again, we consider z as a parameter; consequently, the constant C is independent of z. By applying the 
operator Vjj- to equation (93) we obtain 

-f (R 1 T-0)=0. 

Then, multiplying equation (24) by R\ and adding to the above equation we reach 



(v ■ V H )u -(^J V h - v(x, y, £, t)d£j ^ + (« • V H )v - (V H • v 



(111) 



+Ri 



+ L 1 e + v H ■ 



(v ■ V H )u 



-(£ 



V H • v(x,y,£,t)d£ 



du 



+ (u ■ Vh)v - (Vif • u)u 



92; 



- foV 



1 1 

R2 



(112) 
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Taking the inner product of equation (111) with r/ in L 2 (il) and the equation (112) with 9 in L 2 (fl), integrating 
by parts and observing that r/\ z=a = rj\ z= _ h = 9\ z=a = 9\ z= _ h = 0, we get 



Id 



+ 11*112) , i 



dt 



Ri (llV^Hi + IIV^^Hl) + — {\\d zV f 2 + \\d z 6\\l) 



H 



(v ■ V H )u 



du 



V H ■ v(x, y, £, t)d£ ) — + (u ■ V H )v - (V H ■ v)u 



t] — /o R\Tr] > dxdydz 



} 



L{ WH ( w ' Vff ) u ~(/ VH-v(x,y,t,t)d^^ + {u-V H )v-{V H -v)u 6>J dxdydz 



Ri 



v ■ V H T 



dT 1 



V H -v(x,y,Z,t)dt) _ + - -Q+ -5- -"5- T - 



9z /i 



1 1 



J?3 J?o 



9 \ dxdydz 



< C\\Q\\ 2 \\e\\ 2 + C [ (\v\ \V H u\ + \u\ \V H v\ + \u z \ [ \V H ■ v\ dA (\V HV \ + \V H 9\) dxdydz + C\\T\\, 
Jn \ J-h / 



+C 



\V H v\\T\\6\ + \v\\T\\V H 6\ + \T z \ 



\V H -v\dt) \9\ 



dxdydz + C|M| 2 ||V ff 0|| 2 + 0||T Z || 2 ||^|| 2 



< C||Q|| 2 ||0|| 2 + cj (\v\ \V H u\ + \u\ \V H v\ + \u z \ J° (\6\ + \T\) dt) (\V HV \ + \V H 0\) dxdydz + C\\T\\ 2 \\ V \\ 2 



+C 



\V H v\\T\\9\ + \v\\T\\V H 9\ + \T z \ 



+ \T\)d£)\6\ 



dxdydz + C|M| 2 ||Vh0|| 2 + C||T Z || 2 ||^|| 2 . 



Using Holder inequality, and inequalities (44) and (110), we obtain 

i djm+m + _i_ (l|Vffr?lli + hv^iii) + _l ( ||^||2 + lldM) 

< C (||w||6||V H ti|| 3 + ||w||6||V^||3 + KHJ||V*u*|||||0|||||V*0||J + ||T|U|| Uz || 2 ) (\\V H vh + \\V H 6\\ 2 ) 

+O||T|| 2 |M|2 + O||Vm;||2||T||oo||0||2 + q|«|| 2 ||T|| 00 ||V^|| 2 + 011^11211^21^^112 + ^110011^11211^2 

+o||Q|| 2 ||e|| 2 + c||«|| 2 ||VH^|| 2 + o||r z || 2 ||^|| 2 

< C(\\T\\3 + \\v\\I\\Vhv\\I + ||0||I||Vh0||!) (|| V || a + ||u|| 6 ) (\\V hV \\ 2 + \\V H 6h) + 

+C (|K|| 2 /2 ||V^|| 2 /2 ||tf|| 2 /2 ||V^|| 2 /2 + ||T||oo||^|| 2 + \\T z \\ 2 \\Oh + IMl2||T||oo) (||V^|| 2 + \\V H 6\\ 2 ) 
+0||T|| 2 |H| 2 + C(||V^||2 Halloo + 11^11^11^112) ||f?|| 2 + C||Q|| 2 ||^|| 2 + C||^|| 2 ||V^^|| 2 + C||^|| 2 ||^|| 2 . 
By Young's and the Cauchy-Schwarz inequalities, we have 

d{Ml + mt) + J- (HV^Ili + IIV^HI) + i- + \\d z 9\\l) 

<c{i + \\v\\i + \\u\\t + \\t z \\1 + ||^|||||v^||2) (h\\i + Ml) 

+C\\T\\l + C\\Q\\l + C (Ml + ||u||§ + KH2 + \\v\\l + \\V H v\\l + \\T z \\l) (1 + \\T\\D . 
Thanks to (59), (69), (84), (101), (103), and Gronwall inequality, we have 



l! + l|0|l! + 



/' 



^ (\Whv\\1 + \Nh8\\ 2 2 ) + ^ (\\d zV \\ 2 2 + \\d z 9\\l) 



ds < K 7 , 



where 



K 7 = Ce (K 1+ Kl' 3 +Kl' 3 +Kl)t 



+ \\Q\\i + K'i (K 2 + Kl" + K 1 ^ 3 + K 6 ) 



(113) 
(114) 
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3.8. \\A H v" 2 



H 1 (M) 



||Vjt(V£ • Vz)\\ 2 H i m + || V H (Vh • w z + RiT) \\ 2 mm + ||V ff T||| 1(n) estimates. First, let us 



observe that 

\\7 H v(x,y,z)\ < \\7 H v(x,y)\ 
Therefore, from the above and (107), we have 



\V H v z (x,y,z)\ dz. 



OO < IIV^Hoo + 



,0 

/ \Vhu\ dz 
J-h 



< ||V/rtJ||oo + fli / l|V H /?||oo dz + 

-h 



\\7 H C\ dz 



By applying inequality (40) to V hv and f® h |VjjCI dz we reach 
llVirUlloo < C\\V H v\\m(M) (1 + log + \\V h v\\h^m)) 1/2 > 



/ |V ff CI dz 


< c 


/ |V ff CI dz 


[ 1 + log+ 


[ |V ff CI ^ 


J-h 


OO 


J-h 


m(M) \ 





1/2 



H 2 (M) / 



Applying inequality (41) to Vjj/3, also by (105) and (106), we reach 

llV^Hoo < C (\\V H ■ /3|U + || • /3||oo) (1 + log+ \\V h P\\h*{m)) < CWTW^ (1 + log+ ||T|| H2(M) ) 
Therefore, by (115)-(118), we infer that 

||V^|U<C||V^||H 1 (M)(l + log + ||V^||^(M)) 1/2 + C f [Halloo (1 + log + ||T|| H2(M) )] dz 



+C 



riz 



1 + log + 



H 1 (M) 



J-h 



riz 



1/2 



H 2 {M) / 



< C\\V h v\\hhm) (1 + log+ || V ff tJ|| H2(M) ) 1/2 + C||T|U + log+ ||T|| H2(M) 



/ IVhCI <fo 


f 1 + log+ 


[ IVhCI <fo 


J-h 


H 1 (M) V 


J-h 



1/2 



H 2 (M)j 



< C\\V H v\\ m{M) (1 + log+ ||V ff iJ|| ff2( Af)) 1/2 + C IITI^ (1 + log+ ||A H T|| L2(0) ) 
+C (IhllnMO) + ll<?ll^(n)) [1 + 1°S + (II A^|| L2(0) + ||A ff (9|| i2(n) )] 



1/2 



(115) 

(116) 
(117) 

(118) 



(119) 



3.8.1. || Vij Ajjt;||2 estimates. By taking the Ah to equation (29) and then taking the inner product of equation 
(29) with -A 2 H v in L 2 (M), we reach 

l d\\V H A H v\\ 2 2 1 

Ri 



dt 



2 : -w^ini 



= / A H \{v-Vh)v+[(v-V h )v + {Vh -v)tf + fokxv] ■ A 2 h v 
Jm l j 



dxdy. 
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Integrating by parts and applying (30), we obtain 



ld\\V H A H v 



dt 



+ Ti \\^ H v\\i 



C J M [M |V^| + |V H tJ| |V^tJ|+ J (\v\ \V 3 H v\ + \V H v\ \V 2 h v\) dz + \A H v\j \A 2 h v\ dxdy 
C I \\v\ \\7%v\ + \V H v\ \V 2 H v\ + ( f \v\ dz) ( f (IV^CI + \V\P\) dz 

JM I \J-h / \J-h J 



+ / (iVaCl + |V ff /3|) dz 



-h 



J (|V ff Cl + |V|r/9|) dz) + |A H «|| |A ff tJ| dxdy. 



By applying Holder inequality, (33), (36), 39, (44) and (110) to the above estimate, we obtain 



l d||V g A g U||l 1 
2 dt i?i 



<c(||u||4 ||V^||4 + ||V ff TJ|| 4 ||V ff TJ|| 4 + / \\v\\oo dz [ (|V 3 ff CI + Nh(3\) dz 

I J-h J-h 



+ 



[ (IVhCI + |V*/?|) dz 

J-h 



(\v 2 H c\ + \v 2 H P\)dz 



+ \\A H v\\ 2 



} HA 2 ^| 



<c{\\v\\l\\V H v\\l ||V 3 H TJ||J||Al,TJ||| + ||V ff TJ|| 2 \\V%v\ 



+ \\v\\1\Whv\\So (\\A hV \\ 2 + \\A H e\\ 2 + \\A H T\\ 2 ) 

+ (Wvh + \\e\\ 2 + \\t\\ 2 ) (|| v^|| 2 + 1| v%e\\ 2 + \\v 2 h t\\ 2 ) + \\A H e\\ 2 } \\a 2 h v\\ 2 . 

Thus, by Young's and the Cauchy-Schwarz inequalities, we have 
^M + i-||A^ 
<C(\\v\\ 2 \\V H v\\ 2 + \\V H v\\ 2 ) \\V H A H v\\ 2 

+ (NHIVtfSIU + 11^7111 + \\0\\ 2 2 + \\T\\ 2 ) (\\A HV \\ 2 2 + \\A H 9\\ 2 2 + \\A H T\\ 2 2 ) . 



(120) 
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3.8.2. ||A#T||2 + || VffT z ||2 estimates. By applying the operator Ah to equation (24), and then taking the inner 
product of equation (24) with AhT + T zz in L 2 (Q), we get 

l d(||A g T||| + ||V g r z |||) , 1 
2 



dt 

v ■ V H T 



R 3 

V H -v(x,y,^t)dA (% + \) - Q 



- I A H 
Jn 



AhT dxdydz 



A H v -V H T + 2V H v -V 2 H T 



VhT z dxdydz 

A H (y H -v(x,y,£,t))d£ 



dT 1 

dz h 



AhT dxdydz 



-2 (J ^ V h (Vh ■ v(x, y, £, t))d£j V H T Z - A H Q 
+\7 H u ■ V H T + u ■ \7 H T -\7h(Vh-v)( — + -\-{\7 h -v) V H T Z - V H Q 



VhT z dxdydz 



= ~ J A " v ■ V h T + 2V H v ■ V%T - i ^ A H {V H -v{x,y,Z,t))d^J +A H {V H -v)T 
+2\7 H (V H ■ v) ■ V H T - A H Q] A h T dxdydz 

- f f A H (V H - v(x, y, £, t))d£ T + 2 [ V H (V H ■ v{x, y, £, t))d£ ■ V H T 
Jn U-h J-h 

dT 

Vhv ■ V H T Z + (V h (Vh - v ))-q^+ ^hu ■ V H T + u ■ V 2 H T 



AhT z dxdydz 



-V ff (V H • v) + V H (Vjj ■ u) T + (Vfl- ■ u) V H T - V H Qz 



-j 



V ' hT z dxdydz 



dT 



V H (V ff • v{x, y, £, t))d£ ) + V H (Vh • v) T + {W H ■ v) V H T 



V ' hT zz dxdydz. 
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Thus, 



i«!^Zl±l^ + i(lA^B + |v^.B) 

< C J {\V H v\\V 2 H T\ + \A H v\ \V H T\ + \A H {V H -v)\ + \V 2 H {V H ■ v)\ \T\ + \V H (V H • v)\ \V H T\ 
+ \A H Q\} \A H T\ dxdydz 

\A H (V H ■ v)\ dz\T\ + f \Vh(V h ■ v)\ dz\V H T\) \A H T Z \\ dxdydz 



+C 



-•o 

h 



+C f {\Vhv\\V h T z \ + \Vh(Vh-v)\\T z \ + \Vhu\\V h T\ + \u\ \V 2 h T\ 
Jn 

+ |V ff (Vh • w)| + |Vh(Vh • u)| \T\ + \V H Q Z \} |V ff T z | dxdydz 



+C L \I I |Vff(Vff ' w)l<fe ) \ T z\ + \Vh(Vh-v)\\T\ + \\7 h -v\\V h T\ 



\^hT zz \ dxdydz. 



Thanks to (107)-(109), we obtain 

ld(||A g r|||+||V g T,|||) I 
2 



dt 



+ — (\\A H T Z \\ 2 2 + \\V H T ZZ \\ 2 2 ) 



<C [ (|V g i>||V g T| + [ (\A H (\ + \V H T\)dz\V H T\ + f (\A H 6\ + \A H T\) dz 

Jn I J-h J-h 

+ J } {\V 2 H e\ + \X7 2 H T\)dz \T\ + J° ^\V H 6\ + \V H T\)dz\X7 H T\ + \A H Q\^ \A H T\ dxdydz 



+C 



{\V 2 H e\ + \V 2 H T\) dz\T\ + / (|V g 0| + |V g T|) dz\V H T\ \A H T Z \ \ dxdydz 



+C [ \\\7 H v\\\7„T z \+ f (\V H 0\ + \V H T\)dz\T z \ + \\7„u\\\7„T\ 
Jn I J-h 

+\u\ \V 2 H T\ + J (|V g 0| + |V g T|) dz + (|V g 0| + |V g T|) \T\ + \V H Q Z \} |V ' h T z \ dxdydz 

+C [ \( [ (|V g 0| + |V g T|) dz] \T Z \ + (|V g 0| + \V H T\)\T\ 
Jn l\J-h / 

+ / (\0\ + \T\) dz\V H T\] \V H T ZZ \ dxdydz. 
J-h J 
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Using Holder inequality, and inequalities (44), (45) and (110), we obtain 
ld(||A ff T||!+||V*T z ||!) , l_ niA mll2 , ||V7 m u2 , 



2 J t + -r {^hT z \\ 2 2 + ||VffT„|| 2 , 



< C\\V H vU\A H T\\ 2 + C||A ff CI| 2 /2 ||V ff A H CI| 2 /2 ||T|| ^W^hTWT + C\\T\U\V 2 H T\\l 
+C||V^|| 2 (1 + HTIUJHA^TIIa + C\\V H e\\y 2 \\V 2 H e\\\ ,2 \\T\\ ^ 2 ||V|T||^ /2 + C||A*Q|| 2 ||A*T|| 2 
+C{\\V 2 H e\\ 2 + \\V 2 H T\\ 2 ) \\t\u\a h t z \\ 2 + q|V^||^ /2 ||V^||^ 2 ||T||L /2 ||V|,T||^ 2 ||AHr z || 2 

+c||v^|| 00 ||VHr z ||2 + q|v^|| 2 /2 ||v^||^ 2 ||T z ||^ /2 ||VHT z ||^ /2 

+c||T||^|| V | f r||^ 2 ||r z ||f \\v„t z \\1 /2 + clivals ||VhT|| 6 ||v^|| 2 + chu \\v 2 h t\\ 2 \\v h t z \\ 2 

+C(||V^||2 + ||VffT|| 2 )(l + ||T|| 00 ) 

+C (||V ff ^||^ /2 ||V|0||f + C||T|| ^nv^ni^ 2 ) \\T z \\l /2 \\V H T Z \\\ /2 \\V H T ZZ \\ 2 

+ \\y H Q z h \\VhT z \\ 2 + C(\\V H 0h + IIVffTllaJHTIIoollVffT^lla 
+C\\e\\ 1 2 /2 ||V^|| 2 /2 ||T||V2||v^T||^ 2 ||VhT zz || 2 

< C\\V H v\U\\A H T\\ 2 + \\V H T Z \\ 2 ) + C{\\VHVh + II V ff ^|| 2 ) 1/2 (||A^|| 2 + || A H 6\\ 2 fl 2 ||T||^ 2 || V^T||f 
+C{\\A hV \\ 2 + ||A^|| 2 + ||A ff T|| 2 ) (1 + HTIU ||A ff T|| 2 + C||A ff Q|| 2 ||A ff T|| 2 + C\\\7„Q Z \\ 2 \\\7„T Z \\ 2 

+ (||V^|| 2 /2 ||V^||I + ||V^||!||V 2 ^||I + ||T||i||A ff T||j) ||T||i||A H r||I (||A ff T|| 2 + \\A H T Z \\ 2 ) 
+C7||«||e||V ja -r jB |||||V| l -r JB ||}||^ i rr|H-C7 [llV^lllllVfr^llJ H^llillVirr^lll 

+ (ri|4+||e|| 2 /2 ||V^|| 2 /2 + ||r 2 || 2 /2 ||V ff r z || 2 /2 ) ||T||V2||A H T|| 2 /2 ] ||V*T„|| 2 . 

By the Cauchy-Schwarz and Young's inequalities, we reach 

d(||A g T||l + ||V g r,||l) 1 , 2 2 , 
dt ~Ra " 2 H VhJz2 II2J 

< CIIV^IU (||A„T|| 2 + || V H T Z || 2 ) + C (|| V^|| 2 + || V^|| 2 ) 1/2 (|| A^|| 2 + ||A^|| 2 ) 1/2 ||T|| L /2 ||A ff T||^ 2 
+C(||A ff »7|| 2 + ||A^|| 2 + ||A ff T|| 2 ) HTHoo ||A ff T|| 2 + C\\A H Q\\ 2 \\A H T\\ 2 + C\\V H Q z h \\V H T Z \\ 2 

+ (||V^|| 2 /2 ||A^||} + ||V^||II|A^||! + ||T||i||A H T|||) ||T||i||A ff r||I||A H T|| 2 

+ (HV^IbllA^II, + ||V^|| 2 ||A^|| 2 + ||T||oc!|A ff T|| 2 ) \\T\U\A H T\\ 2 + (||V ff T z || 2 + ||A ff T|| 2 ) 

+C [||V^|| 2 ||A^|| 2 ||T,|| 2 ||V ff r z || 2 + (||T|| 2 + ||0|| 2 ||V*0|| 2 + ||i;|| 2 ||V ff T z || 2 ) ||T|U||A ff T|| 2 ] 

< C||A ff Q|| 2 + C||V ff Q z || 2 + C\\T\t + CUV^IU (||A ff T|| 2 + ||V ff T z || 2 ) 

+C (1 + \\VhvW 2 + \\V H e\\ 2 + \\T\\l + \\u\\l + \\T Z \\ 2 ) (\\A H T\\ 2 + \\V H T Z \\ 2 + \\A hV \\ 2 + \\A H 0\\ 2 ) . 
Thus, we get 

rf (||A g T||l + ||V g T 2 ||l ) + _^ (||A ^ ||2 + ||VA||2) 

< C||A H Q|| 2 + C||V ff Q z || 2 + C\\T\\t + C||V^|U (||A ff T|| 2 + ||V ff T z || 2 ) 
+C (1 + IIV^H 2 + ||V^|| 2 + \\T\H + \\u\\l + \\T Z \\ 2 ) (\\A H T\\ 2 + \\V H T Z \\ 2 + \\A HV \\ 2 + \\A H 6\\ 2 ) . (121) 

3.8.3. ||Vij(V# -v z )\\ 2 H1 ^ + \\Vh (Vjj • v z + RiT) \\ 2 H i^ estimates. By acting with A H on equation (111) and 
equation (112), then taking the inner product of equation (111) with Anr] + r\ zz in L 2 and equation (112) with 
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A H 9 + 9 Z 2 in L 2 , respectively, we get 

1 d (HAgrylll + ||V g ^||l + ||A g g||l + ||V g <y 1) 

2 dt 

+i- (HVhA^HI + \\VhAhVz\\1 + ||V ff A H 9\\ 2 2 + \\A H 9 Z \\ 2 2 ) 



+±- (\\A HVz \\ 2 2 + \\VhVzz\\1 + \\A H Z \\ 2 2 + HVh^IH) 



= v H <v 



H ■ 



{v ■ V H )u - Q ^ V„ • v(x, y, C t)d^j ^ 



+(« • V ff )w - (Vh • v)u] + f T - 0)} • V H (A H V + Vzz) dxdydz 

+ L Vh { Vh ' (w ' Vh)u ~{J , Vh ' y ' ^ £ 

+(« • V H )v - (V H • v)u] + f v} • V H (A H + 9 ZZ ) dxdydz 
- jiZx Vif« • V H T + v ■ V 2 H T - (y V H {V H - v)(x, y, £, t)d^j ( 



<9t i 

9z /i 



X _C/- 



1 1 

R3 R2 



V H T Z . 



V H (A H + 6 zz ) dxdydz 



V H (V H ■ v)(x, y, C t)d£ V H T Z A H 9 + 



[ V H -v(x,y,Z,t)dZ)V H T z -V H e z ) 

J-h 



dxdydz 



C [ {(\u\ \V H v\ + \V H u\ \V 2 H v\ + \V 2 H u\ \V H v\ + \V H u\\v\ + \u z \ [ \V 2 H (V H • v)\ dz 

JO I V J-h 

+ \V H u z \ J \V H (V H ■ v)\ dz + \V 2 h u z \ J \V H -v\dz^ (\VHAHv\ + \VHVzz\ + \VHA H e\ + \V H e zz \) 
|Virt;||VffT| + |w||V|rr| + (1 + f \V H {V H ■ v)\ dz + \V H T Z \ f \V H ■ v\ dz 

J-h J-h 

ik'lk |VAl (l V » A **l + l V **"l)} dxdydz + C{\\A H T\\ 2 2 + \\A H v\\ 2 2 + \\A H 0\\ 2 2 ) 
f (|V 3 ff CI + IVf^l) dz + (|VhCI + |V ff /3|) f {\V 2 H C\ + \V 2 H P\) dz + (IV^CI + \V H P\)H 

J-h J-h 

+(1^(1 + 1^/31) f (\V H (\ + \VH(3\)dz + \u z \ [ {\V 2 H e\ + \V 2 H T\) dz 

J-h J-h 

+(|VtfCzl + |Vff/y) / {\v H e\ + \v H T\)dz 

J-h 



+c 

+ \Vr(}\ r 

< c 



(\v H A H v\ + \Vhvzz\ + \v H A H e\ + \v H e zz \) 



+(|V^C*I + |V^|) / (\0\ + \T\)dz 

J-h 

+C \V H T\ / (\V H (\ + \V H (3\)dz+\v\\V 2 H T\ 

J-h 

+(1 + \T Z \) f (\V H 0\ + \V H T\)dz + \V H T z \ [ (\6\ + \T\)dz 

J-h J-h 



+|V ff Q| 



R3 R2 



\V H T Z 



(\V h A h 8\ + \Vh0zz\)\ dxdydz + C(\\A h T\\ 2 + \\A hV \\ 2 + \\A h 9\\ 2 ) . 
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Using Holder inequality, and inequalities (44), (45) and (110), we obtain 

1 d(\\A H r,\\ 2 2 + HVgTfclH + ||A g fl||l + HVgflslH) 

2 dt 

+ -L (\\V H A HV \\l + \\V H A H r, z \\ 2 + \\V H A H 6\\ 2 + \\A H 6 Z \\ 2 ) 

Hi 



+4-(l|Agr? 

-fl2 



z\\2 ■ 



< C(\\V H A HV \\ 2 

+ 



W^HVzzh 



\\A H d z \\ 2 + \\V H 9 zz \\ 2 ) 

\\v H A H e\\ 2 + \\v H e zz \\ 2 ) [(\\u\u + \\v\u (nv^il 



\V 2 H 6\\ 2 



V 2 ff T|| 2 ) 

(NIIIuv^II! + ||f||II|v^||j + ||T|u) (liv^lllliv^lll + liv^lllliv^ni + ||r||i||v 2 ff r||i) 

+ (lIV^IllllV^HI + IIV^IllllV^HI + ||T||i||V 2 H T|||) (hllJllV^IIJ + ||0||!||V ff 0||| + ||T|U) 
+ (ll^lllllV^II! + millllV^Hl) (IIV^IIIIIV^HI + ||T||i||V ff T||l) 
+\K\\I\\Vhu z \\I (lIV^IIIIIV^III + ||T||i||V ff T|||' 



+c 
+ 
+ 
+ 



(liv^lllliv^iil + llv^Hiiiv^iij) (ll^lllliv^ll! + ||T|U) 

(IMIIIIV^H! + lieillUV^IIJ + ||T|| 4 ) ||T||i||AgT||I + \\v\U\A H T\\, 
(IIV^IIIIIV^III + ||T||i||A ff r|||) \\T Z \\I\\V H T Z \\I 

(» 



T|U \\V H T z \\i\\A H T z 



llliv^ll! 

R3 R2 

By Young's inequality and Cauchy-Schwarz inequality, we have 



Vu0||2 + ||V ff T|| 2 + ||VjrQ||2 



zz\\2 



(\\\7HA H e\\ 2 + \\v H e zz h) + c(\\A H T\\l+ \\a hV \\ 2 + \\A H e\\ 2 ) 



d {WAnnW 2 + \\VhVzWI + \\A H e\\ 2 + \\v H e z \\ 2 ) 



dt 



ft 



(\\VhA hV \\1 + \\VhA hVz \\ 2 2 + \\V H A H 9\\ 2 + \\A H 6 Z \\ 2 ) 



i?2 



(IIA^Hi + HV^ z 



Ah^zIIo 



\V H 0z 



< C (\\A H T\\ 2 + \\V H T Z \\ 2 + \\A HV \\ 2 + \\VhVz\\1 



IIA^II 2 + ||V ff ^|| 2 ) [1 + \\T\\t + \\\7 HUz \\ 2 + \\u zz \\ 2 



-\\Vhv\\ 

2 



ll^lli 



\m\l 



liv^ll 

+ hll2l|V^||^ + ||0||2||V^|| 2 + ||T|| 4 o 
R\{R\ + J?2)(-R2 — -R3) 2 



\T, 



\Vhv\\ 



\\Vh6\\1] 



IKIII'llV^III' + llVtfQII 2 



+ 



E>2 I?2 
K 2 H 3 



-\\VhT z . 



(122) 



Next, we will obtain an estimate for 

\\V H A H v\\ 2 2 + ||V«(V H • v z )\\ 2 Hl{n) + || Vg (Vh • v z + KT) \\ 2 Hl(n) + C R \\ V ff T|| 2 ffl(0) , 
where Cr = ^(fii+gX^-fla) 2 _ Dcnote by 

X = 1 + ||Vh A H v\\ 2 2 + C R \\A H T\\ 2 2 + C R \\V H T Z \\ 2 2 + \\A HV \\ 2 2 + ||V ff 7? z || 2 + \\A H 9\\ 2 2 + \\V H 0z\\ 2 2 , 
y = \\A 2 H vf 2 + \\A H T z f 2 + \\V H T ZZ \\ 2 2 + \\V H A hV \\ 2 2 + \\A HVz \\ 2 2 + \\V H A H 9\\ 2 2 + \\A H 6 Z \\ 2 2 . 
Thus, by (119), we get 

IIV^Hoo < C {\\S7 H v\\ HHM) + WT^ + \\v\\ H ^n) + Mh^)) log A". 
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By virtue of (120), (121), (122), Young's inequality, Cauchy-Schwarz inequality and the above, we obtain 
dX 

_ + C y < c\\v\\ 2 (||v ff TJ|| ffl(M ) + + |M|„ 1(n) + \\9\\ H1(n) ) x \o g x 



+ 



\L+ \\T z g + \\v\\ 2 2 (l + \\v\\h l ,r 

|2\ |i II 2 i It i ||/)||2\ i|/)||2 



+\K\\%, + (i + IMIi)ll* + (i + IIO PWIA x 

+C{\\A H Q\\l + \\V„Q z \\l + ||r,||i||V ff »7||i + H^llHV^ + \\T\t + H^lllllV^Hi} . 
Let X = e z . Then, f = e 2 f . As a result we have 

<IZ <C\\v\\ 2 (l|VHtJ||H 1 (M) + im| 00 + ||r ? ||H 1 (a) + ||^||^ ( o)) -Z 



eft 

+ 



1 + \\A H Q\\l + \\V H Q Z \\ 2 2 + \\TWl + \\T z \\l + \\v\\l (l + \\v\\ 2 HHMh 

+ (1 + \\u*\\l) \\u z \\h + (1 + ll^lli) \\v\\ 2 m + (1 + II* II!) Il*l&i] • (123) 
Thanks to Gronwall inequality, and the estimates established in the previous subsections, we get 

Z<K, X <e K (124) 

where 

K = e c(K 1+ K 2+ K 7 ) ^ + || Uo|| ^ 4 + j|To||2ff2 + t + \\ Ah q\\1 t + \\v H Q z \\l t] . (125) 
Moreover, we have 

rt 

yds<K s , (126) 



/ 

Jo 



10 

where 

Ks = e c(K 1+ K 2+ K 7) ^ + n^u^ + ||To||2ff2 + t + \\ AllQ f 2 t + ||V H Q Z ||| t] . (127) 

Thanks to (103) and the above we conclude that the quantities J Q * ||^zz(s)||#i(Q) ds, J Q * ||AjjVjjt) z (s)||| r i^ rfs, 
II A/fW^(t)|| ff i( f2 ), ||VijT(t)|| H i(n), and J Q || VhT z (s)\\ 2 h1 ^ ds are all bounded uniformly in time, t, over the interval 
[0,%). Therefore, the strong solution (v(t),T(t)) exists globally in time. □ 

4. Uniqueness of the Solutions 

In this section we state and prove the global existence and uniqueness of the strong solution of system (22)-(28). 

Theorem 2. Suppose that Q e H 2 (Q). Then for every v G H 4 (Q), T e H 2 (Q) and T > 0, there is a unique 
solution (v,p s ,T) of system (22)-(28) with 

A H v z , V H T eC([0,7],-ff 1 (ft)) ) 

v zz , A H \7 H v z , W H T Z ^L 2 {[Q 1 T],H 1 (^)), 

Proof. In Theorem 2 of the previous section we proved that the strong solutions exist globally in time. Therefore, 
it remains to prove the uniqueness of strong solutions, and their continuous dependence on initial data, in the sense 
specified by equation (140) below. Let (v\, (p s )i,Ti) and (v 2 , (p s ) 2 ,T 2 ) be two strong solutions of system (22)- 
(28) with initial data ((«o)i, (7o)i) and {{v ) 2 , (T ) 2 ), respectively. Denote by <j> = Vi -v 2 ,q s = {p s ) x - {p s ) 2 ,ip = 
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T\ — T 2 . It is clear that 

d<t> 



dt 



Licfi + (vi ■ Vh)4> + {<$> ■ Vij)«2 - / V H -Vi&y,^^ 



dz 



<i>{x,y,£,,t)d£, 



dv 2 
dz 



+f k x 4> + V H q s - V H i/>(x, y, £, t)<%j = 0, 



at it 3 



(/, ■ v 1 (x,y,U)d^j ^ - ■ (t>(x,y,^,t)d^j (j^ + ^J =0, 



with initial data 



<p(x,y,z,0) = (v )i - (v ) 2 , 
iP(x,y,z,0) = (T ) 1 -(T ) 2 . 
Taking the inner product of equation (128) with <fi in L 2 (Q), and equation (129) with tp in L 2 (Sl), we get 



(128) 



(129) 

(130) 
(131) 



1 



2 dt i?i 



: 1 ^iiv^ + -^,i-; 



= -/ 



(vi • Vh)<P + {<i> ■ Vh )V2 - ( / V ff -vi(a;,i/,^t)de 



+ 
and 



fok x <j> + V H q s -V H [ / i>{x,y,£,i)d€ 

h 



dz 
dxdydz, 



"(/> 



<fi(x,y,£,t)d£ 



dv 2 
~8z~ 



dxdydz 



1« , 1 



ll^lll 



2 rfi i? 3 

V ff • vi(x,y,£,t)d£ 



[v! ■ V H V + • V H T 2 



5^ 



\dz h 



ip dxdydz. 



Notice that 

f Q kX(f>-tp = 0. 

Integrating by parts, and using the boundary conditions (25) and (26), we have 



/ n (e>,.v„) 



V// ■ /•.!.(. r. //.(./ >< ) ^ ) ■ <> d.rdydr. = 



~X ( Wl ' Vir ^~ (/ Vir-t;i(a;,|/,^t)^ dxdydz = Q. 

Integrating by parts, and using the boundary conditions (25) and (26), we get 



(132) 

(133) 
(134) 



dxdydz 



f ( f il>{x,y,£,t)dn (V H -(p) dxdydz. 

JQ \J-h J 



(135) 
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Thus, by (132)-(135) wc have 
1« 



2 dt 



ri\ il 2 



(</> ■ V h)v2 ■ 4> dxdydz 



and 



J J V ff • cj>(x,y,£,t)d£^ ■ <p dxdydz + J (^J ip(x,y,^,t)d^j (V H ■ 4>) dxdydz. 

\d\m 



2 dt 



1 

R~3 



IIV'.Ill 



VhT 2 )iP dxdydz 



- [ (Vjj • </>) T 2 ?/> dxdydz - / / Vj? ■ y, £, t)d^T 2 ^ z dxdydz. 
Jvi Jn J-h 

Notice that by Holder inequality and (44) 



• V# )«2 • 4> dxdydz 
VifT 2 V dxdydz 



< 



M |0| |V.ff0| dxdydz 



< Ihllell^llJIIV^H 



3/2 
2 : 



y y Vh ■ <j)(x 7 y,£ 7 t)d£^- ■ < ■ dxdydz 



< ||VhT 2 || 3 ||V||2||0||6 < C||VhT 2 || h1 ||V|| 2 ||V^|| 2 ; 

<% 2 



<C||V ff 0|| 2 



a« 2 



9z 



ll*<c 



dv 2 



dz 



12 i 



/ I* V H -<j>(x,y, UH^ip dxdydz < [ [ \V H cj>\dz [ 

JnJ-h 0z JO. J-h J-h 

f f \V H <t>\dz ( f 

JnJ-h \J-h 



dT 2 



dz dxdy 



< C 



< c 



dzl \ip\ dz) dxdy 

||V^||2||V||2. 





dT 2 


2 


/ 




dz 


J-h 


dz 



1/2 



l|Vi^|h 



< c 



9AhT 2 



dz 



Therefore, by estimates (136)-(139), we reach 

< C (||«2|| 6 ||0IIII|V^||2 /2 + \\V H T 2 \\ H1 W\\ 2 \\VHCPh + ll^2|| 6 ||0||!l|V^||2 /2 )+C 

By Young's inequality and the Cauchy-Schwarz inequality, we get 



dA H T 2 



dz 



(136) 
(137) 

(138) 



(139) 



liv^lh 



Thanks to Gronwall inequality, we obtain 

\\m\l+\\m\\l<(\m=o)\\i+\mt=o)\\i)x 



8A H T 2 



dz 



(IHIl + Wl). 



cxp < C 



\v 2 (s)\\t+\\V H T 2 (s)\\ 2 H1 + \\d z v 2 (s)\\i + 



OA H T 2 



dz 



ds 
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As a result of (124), we have 

\W)\\l + \\m\\l < (\W = 0)||1 + HVC* - 0)111) cxp [C ((K 2 3 /3 + K + KT) t + Ks)}. (140) 

The above inequality proves the continuous dependence of the solutions on the initial data, and in particular, 
when <j){t = 0) = and ip(t — 0) = 0, we have <f>(t) = and ip(t) — 0, for all t > 0. Therefore, the strong solution 
is unique. 
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